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Abstract 



We construct quasimodes for the Klein-Gordon equation on the black 
hole exterior of Kerr-Anti-de Sitter (Kerr-AdS) spacetimes. Such quasi- 
modes are associated with time-periodic approximate solutions of the 
Klein Gordon equation and provide natural candidates to probe the decay 
of solutions on these backgrounds. They are constructed as the solutions 
of a semi-classical non-linear eigenvalue problem arising after separation 
O ' of variables, with the (inverse of the) angular momentum playing the role 

of the semi-classical parameter. Our construction results in exponentially 
j_l ■ small errors in the semi-classical parameter. This implies that general 

solutions to the Klein Gordon equation on Kerr-AdS cannot decay faster 
than logarithmically. The latter result completes previous work by the au- 
thors, where a logarithmic decay rate was established as an upper bound. 
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1 Introduction 

There is currently a lot of mathematical activity concerning the analysis of waves 
on the exterior of black hole backgrounds!]] The main motivation is the black 
hole stability problem, i.e. the conjectured non-linear asymptotic stability of 
the two-parameter family of asymptotically-flat Kerr spacetimes (M, <?M,a), the 
latter being stationary solutions of the vacuum Einstein equations Ric [g] = 0. 
From the point of view of non-linear partial differential equations, the analysis of 
linear scalar waves on black holes is a prerequisite to successfully understand the 
non-linear hyperbolic Einstein equations in a neighborhood of the Kerr family. 

A major insight that crystallized in the last decade [T2l H4l [T51 [TBI I3"2l 1371 fel 
[T71HB] is that the fundamental geometric obstacles to the decay of waves, namely 
superradiance and trapped null-gcodcsics, can be overcome by exploiting the 
normal hypcrbolicity of the trapping, the redshift effect near the event horizon 
and the natural dispersion of waves in asymptotically flat (and asymptotically dc 
Sitter) spacetimes. In particular, polynomial decay rates have been established 
for solutions to the scalar wave equation on the exterior of any member of the 
sub-extremal Kerr family of spacetimes [T5] . 

Changing the black hole geometry can have dramatic effects on the behav- 
ior of linear waves through the subtle interplay of the redshift, the superradi- 
ance and the trapping. Aretakis [3J 0] showed that for extremal black holes 
(whose vanishing surface gravity leads to a degeneration of the redshift effect) 
the transversal derivatives of general solutions to the wave equation will grow 
along the event horizon. In |35| . it is proven that for the massive wave equation 
on a sub-extremal (|a| < M) Kerr spacetime, exponentially growing solutions 
can be constructed on the exterior, exploiting an amplification of the superra- 
diance caused by the confining properties of the mass term. 

In this paper, we shall be interested in a black hole geometry for which a 
strong trapping phenomenon leads to a very slow (only logarithmic) decay of 
waves. More precisely, we will study the behavior of solutions to the massive 
wave equation 




(1) 



1 See 1131 for an introduction and a review of recent results. 
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in the exterior of asymptotically Anti-de-Sitter (AdS) black holes with spacetime 
metric g. 

Due to their AdS asymptotics, these spacetimes are not globally hyperbolic. 
Nonetheless, the equation (P) is well posed in suitably weighted Sobolev spaces, 
denoted here by H\ dSl provided a satisfies the Breitenlohner-Friedmann bound 
a < 9/4. See [25] as well as [351 H3 G3 125 and [SS] for a complete treatment of 
general boundary conditions. 

The global properties of solutions to (H} on the exterior of non-superradiant@ 
Kerr- AdS black holes were studied in [24 , 26, 30 . In particular, boundedness was 
obtained in [24[ [30] and logarithmic decay in time for general H\ dS solutions in 
|26j . We summarize these results in the following theorem. We refer to Section 
I2.1.1I for the precice definitions of the Kerr- AdS spacetimes, the area-radius r + 
of the event horizon and the S t » -foliation and to Section |2~21 for the definitions 
of the norms and energies used in the statement below. At this point we only 
remark that e\ [ip] is an energy density involving all first derivatives of ip while 
e 2 [ip] and e 2 [tp] involve all second derivatives (with appropriate weights): 

Theorem 1.1. Let (g,TZ) denote the black hole exterior of a Kerr-AdS space- 
time with mass M > 0, angular momentum per unit mass a and cosmological 
constant A = — 75. Assume that the parameters satisfy a < j, \a\ < I. Fix a 
spacelike slice S t * intersecting TL + . Then the following is true. 

1. Equation rijj is well-posed in CH\ ds on (g,7Z) for any k > 2 for initial 
data prescribed on S t * . See I25f . 

2. The solutions of fl]) arising from data prescribed on E t * remain uniformly 
bounded on the black hole exterior provided r 2 ^ > \a\l holds. In particular, 

I ei [ip] (t*) r 2 sin OdrdOdc/) < [ e x [ip] (t* ) r 2 sin OdrdOdcP . (2) 

Analogous statements hold for all higher H^ dS -norms. In particular, 

I e 2 [ip](t*) r 2 sin ddrdOdcp < j e 2 [ip] (*q) r 2 sin 6drd6d<P . (3) 

and the same statement for e 2 [ip] (t*). See H^H 

3. If the parameters satisfy r\ > \a\l, the solutions of $7$) satisfy the following 
global decay estimate: 

ei [ip] (t*) r 2 sin BdrdBd<p < ^ / e 2 bP] (*o) 7,2 sm6drd6d<p 

pog(2 + t*)] Js t j 

(4) 

for all t*>t*>0. See 

2 This means that the parameters of the black hole satisfy > \a\l. See Remark 1 1 . 2 1 and 
Section [2] 

3 The aforementioned papers (as well as 1251 ) are only concerned with the e2 [i/ , ]-energy. It is 
remarked that by commutation with angular momentum operators one can prove boundedness 
for the e2 [?/>]-energy (which differs from the (=2 [?/<]-cncrgy through the weights of the angular 
derivatives). For completeness, we provide an explicit proof of this statement in the appendix 
of this paper. 
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Remark 1.1. The constant implicit in the symbol appearing in\^ and\3\ 
depends only on the fixed parameters M , i, a and a. 

Remark 1.2. The condition r\_ > \a\l on the parameters in\^ and\^ guaran- 
tees the existence of a globally causal Killing vectorfield on the black hole exte- 
rior, the Hawking- Reall vectorfield }23\/ , which explains why such black holes are 
sometimes referred to as "non-superradiant" . If one restricts to axisymmetric 
solutions of (QP, this condition can be dropped for both and [3J in Theorem 

o 

Remark 1.3. Note that the boundedness statement fD) does not lose derivatives, 
while the decay statement ^ does. This is the familiar loss of derivatives caused 
by the existence of trapped null- geodesies \3J$ . 

This logarithmic decay rate Q was conjectured to be sharp in (53] in view 
of the discovery of a new stable trapping phenomenon, itself a consequence of 
the coupling between the lack of dispersion at the asymptotic end and the usual 
(unstable) trapping on black hole exteriors. See again [25] , 



which is determined in Lemma 13.11 as the location of the unique 



I. 1 The main results 

In this paper, we shall prove that the logarithmic decay estimate of Theorem 

II. H is indeed sharp. Recall that for the obstacle problem, it is classical [Mj that 
lower bounds on the rate of energy decay can be obtained from the construction 
of approximate eigenfunctions, also called quasimodes, of the associated elliptic 
operator, obtained by formally taking the Fourier transform in time of the wave 
operator. Our main theorem establishes the existence of such quasimodes for 
Kerr-AdS with exponentially small errorsQ 

The statement of the following theorem will involve the quantity r max S 

3M 

maximum of a simple radial function. For Schwarzschild-AdS, r max = 3M will 
be the location of the well-known photon sphere. 

Theorem 1.2 (Quasimodes for Kerr-AdS). Let (g,7V) denote the black hole ex- 
terior of a Kerr-AdS spacetime, with mass M > ; angular momentum per unit 
mass a and cosmological constant A = — Jy. Assume that the parameters satisfy 
a < j, \a\ < I. Let (t,r,0,<p) denote standard Boyer-Lindquist coordinates on 
1Z. Then, for S > sufficiently small (depending only on the parameters I, M, 
a, a ), there exists a family of non-zero functions %pe £ H\ ds for any k > such 
that 

1. ipz(t,r,6,ip) = e lulet (p£(r,9) (axisymmetric and time-periodic), 

2 

2. < c < gre+i) ^ C > f or constants c and C independent of £ 
( uniform bounds on the frequencies ), 

3. for all t* > t*, for all k >0,\\ (D g + £) M H * ds (s t *) < C k e- c « l \\^\\ H%s[ ^ Q) , 
for some Ck > independent of I (approximate solutions to the wave equa- 
tion), 



4 Note that in order to deduce the sharpness of the logarithmic decay rate from the quasi- 
modes, polynomial errors would a priori not be sufficient. 
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4- the support of Fg := (O g + £r) tp£ is contained in {r max < r < r max + 6} 
(spatial localization of the error), 

5. the support ofipe(r,6) is contained in {r > r max } 
(spatial localization of the solution). 

Note that the ipg have constant -ff^ dS -norms and hence exhibit no decay. On 
the other hand, a standard application of Duhamel's formula shows that the tpt 
are good approximations to the solution of (O g + j?) "ip = arising from the 
data induced by ipt, at least up to a time t ~ e h . 

Corollary 1.1. Let (TZ,g) denote the black hole exterior of a Kerr-AdS space- 
time as in Theorem ] 1.21 Denote by SCH\ dS the set of CH\ dS solutions to ([T]) 
with a < j. Let £q > be fixed and define for any non-zero ip and t* > 



W (**) :=log(2 + i*) 



/E t .n{r>r ma:r } 



d [tp] (t*) r 2 sin 6drd6d(f> 



f s t e 2 [if>] (^) r 2 sin 6drd6d(f> 



Then there exists a universal (depending only on M, a, \a\ and I) constant 
C > such that 

limsup sup Q [ip] (i*) > C > 0. 

Remark 1.4. Corollary [T7T| implies that the (semi-)local energy in {r > r max } 
cannot decay universally faster than (logt*) , unless one loses more deriva- 
tives. 

Remark 1.5. We emphasize that no smallness assumption on the angular mo- 
mentum a is needed apart from the condition \a\ < I which ensures that the 
metric is a regular black hole metric. 

Remark 1.6. Note that the ipt constructed in Theorem ] 1.2\ are axisymmetric, 
while the decay estimate of Theorem \1.1\ holds for general solutions (provided 
r^_ > \a\l holds in the non- axisymmetric case). Since we are concerned here 
with a lower bound on the uniform decay rate, an analysis within axisymmetry 
is sufficient. This allows us to drop the non- super radiant condition > \a\l 
in the analysis. On the other hand, for (sufficiently) superradiant black holes 
r+ < \a\l one can adapt the proof of \35f to construct exponentially growing 
solutions. Hence in this case, the quasimodes we construct are not the "worst" 
solutions on these backgrounds. 

Remark 1.7. Note that the stable trapping occurs only in the region r > r max 
and is associated to certain frequencies. As a consequence, stronger local energy 
decay in r < r max or for some frequency projections of solutions is a priori 
compatible with the results of this paper. 



1.2 Related works and discussion 

1.2.1 Non-linear analysis on asymptotically AdS spacetimes 

In [27l [29], the non-linear spherically-symmetric Einstein-Klein-Gordon system 
for asymptotically AdS initial data was studied, and in particular, the asymp- 
totic stability of Schwarzschild-AdS was proven within this model. 
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For a discussion connecting the logarithmic decay to the non-linear stability 
or instability of asymptotically AdS black holes, we refer to Section 1.4 of [2"rj] . 
We also mention the recent heuristic analysis of |33| drawing attention to a po- 
tential stability mechanism caused by the lack of exact non-linear resonances 
in this setting. For AdS itself, instability was conjectured in [TJJ E] an d [10] - 
More recently, both numerical and additional heuristic evidence has been pre- 
sented |7j . Finally, let us note that asymptotically AdS solutions to the Einstein 
equations have been constructed in |21j . 



1.2.2 Quasi-normal modes of the asymptotically AdS black holes 

Quasi-normal modes, also called resonances, are complex frequencies general- 
izing the well-known normal modes to systems which dissipate energy. There 
is a strong connection between quasimodes and resonances |36j . One way to 
mathematically define them is as poles of the meromophic continuation of a 
truncated resolvent. In the case of asymptotically de- Sitter black holes, this 
theory has been very successfully developed, see [8] [17j [18] . In a recent paper 
[22], Gannot has established, in the case of Schwarzschild-AdS, the existence 
of a sequence of quasi-normal modes (based on an independent construction of 
quasimodes) , indexed by angular momentum £, and with imaginary parts of size 
O (exp(— ^/C)). In particular, his construction confirms the numerical results of 
[T5] and provides an independent proof of Theorem 11.21 and Corollary 11.11 albeit 
restricted to the Schwarzschild-AdS case. While we restrict ourselves here to 
the construction of quasimodes, we strongly believe that our results can be used 
as a basis for the construction of resonances in the Kerr- AdS case for the whole 
range of parameters satisfying \a\ < I and r\ > \a\l. 

1.2.3 Universal minimal decay rates of waves outside stationary black 
holes 

In the context of the obstacle problem in Minkowski space, a celebrated result 
of Burq [§] establishes a logarithmic decay rate for the local energy of waves, 
independently of the geometry of the obstacle causing the trapping. For waves 
outside black holes, in view of the results of [2S], a natural conjecture is: Given 
any black hole exterior (TZ, g) of a stationary spacetime, a logarithmic decay 
of energy similar to that of [26] will hold, provided a uniform boundedness 
statement is true for solutions to the wave equation Dgip = on 1Z. 



1.3 Outline and overview of the proof 

Section [21 introduces the family of Kerr- AdS spacetimes as well as the norms re- 
quired to state our estimates. In Section [3] we exploit the classical fact that the 
wave equation separates on Kerr- AdS. An ingredient which considerably simpli- 
fies our analysis here is the important observation that we can restrict ourselves 
to axisymmetric solutions. With axisymmetry, the separation of variables leads 
to relatively simple, one-dimensional, second order ordinary differential equa- 
tions for the radial functions. In the case of Schwarzschild, they are roughly of 
the form of the semi-classical problem 



G 



for a potential V a (r) , whose general form is depicted below[f| In Section 13.31 
we shall describe in detail the analytic properties of the potentials appearing in 
these equations. 




To construct the quasimodes, we first construct eigcnfunctions for the problem 
§5§ with Dirichlct conditions u = imposed on u at r = r max and r — > oo 
(Section [4]). In particular, we prove a version of Weyl's law, ensuring that for 
any energy between p < E < V max = Va (r max ) we can find (lots of) eigenvalues 

gp^n of ||SJ) in a strip [E — 5, E + 5] . In the Kerr case, the eigenvalue problem 
§5§ turns into a problem of the formal 

1 uj 2 

which is non-linear in lo 2 . An application of the implicit function theorem 
together with global estimates on the behavior of the eigenvalues still allows 
us to conclude the existence of eigenfunctions up of ([S]) with corresponding 
eigenvalues in the range (jj, E + <5] . These estimates, together with the analysis 
of the potential in the Kerr-AdS case, constitute the core of our paper. 

In Section [SJ we recall the so-called Agmon estimates, their proof being 
included to make the paper self-contained. These estimates quantify that the 
solutions ue constructed from the above eigenvalue problem decay exponentially 
in I in a region [r max , r max + 5'}. 

In Section [6l the quasimodes are constructed by cutting off the solution ue 
of the eigenvalue problem in [r max , r max + 5'] so that it vanishes with all deriva- 
tives at r max and then continuing it to be identically zero in [rh OZl r max 1. The 
function <f>e thus constructed will be defined on [r^oz, oo) and the corresponding 
wave function ipe = e lLJe 4'£ wm satisfy the wave equation everywhere except in 
the small strip [r max , r max + 8'] , where the error is exponentially small by the 
Agmon estimate. 

5 For the purpose of this exposition, we neglect terms of lower order in jrfpn m * ne 
potential, as well as the mass-term. The latter is actually unbounded and needs to be absorbed 
with a Hardy inequality. We suppress such technical difficulties in the present discussion. 

6 The fie (a 2 u 2 ) generalize the familiar spherical eigenvalues 1(1+1) of the Schwarzschild 
case to Kerr. See Section |3. II 
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In the last section, we prove Corollary 1 1 . 1 1 using the Duhamel formula. Fi- 
nally, Appendix [A] contains a proof of boundedness for the second energy used in 
this paper. This boundedness statement differs from that obtained in [531130] in 
that it allows for stronger radial weights near infinity for the angular derivatives. 
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2 Preliminaries 

2.1 The Kerr-AdS family of spacetimes 

We recall here some basic facts about the family of Kerr-AdS spacetimes re- 
quired in the paper. We refer the reader to the detailed discussion in our [26] . 

2.1.1 The fixed manifold with boundary 1Z 

Let 1Z denote the manifold with boundary 

K = [0,oo)xlx§ 2 . 

We define standard coordinates y* for [0, oo), t* for R and (6, 4>) for S 2 . This 
defines a coordinate system on 7Z, which is global up to the well-known degen- 
eration of the spherical coordinates. We define the event horizon T-L + to be the 
boundary of 1Z: 

U + =d1Z = {y* = 0}. 

The manifold 1Z will coincide with the domain of outer communication of the 
black hole spacetimes including the future event horizon T-L + . 

2.1.2 The parameter space and the radial function 

A Kerr-AdS spacetime is characterized by its mass M > 0, its angular momen- 
tum per unit mass a and the value of the cosmological constant A < 0. For 
convenience, we shall use mostly I = y^fp instead of A. For the spacetime to 

be regular, we will require \a\ < I. Let us thus fix M, I > and \a\ < I. We then 
define r + (M, I, a) > as the unique real solution of A_(x) = where 

A_{x) = (^1 + ^ (x 2 + a 2 )~2Mx. (7) 

We now define a function r on 1Z as follows. As a function of (t* ,y* ,(p,9), 
r only depends on y* and is a diffcomorphism from [0,oo) to [r + ,oo). The 
collection (t*, r(y*), cp, 9) then form a coordinate system on 1Z, global up to the 
degeneration of the spherical coordinates. Moreover, the horizon H + coincides 
with {r = ?'+}. 
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2.1.3 More coordinates: r*, t, <f> 

Let us define r* by 

dr* , , r 2 + a 2 , 

where A_(r) is given by Q. Note that r*(r + ) = — oo. 

By a small abuse of notation, we shall often write for functions / and g, 
f(r*) = 9{r), instead of f(r*) = g{r(r*)) or f{r*{r)) = g(r). 

Finally, let r cut = r + + 1ma ' '+ , with r rnax > r + defined in Lemma 13.11 
depending only on the parameters M, a and I, and x( r ) be a smooth cut-off 
function with the following property 



1 if r G 



r+, r -^^ 



[0 if r > rcut. 

We introduce the time coordinate t and another angular coordinate <j) as 
t = t* — A (r) x ( r ) an d = </> — B (r) x (r) , 

where 

cL4 _ 2Mr dB _ a (l ~ 7*") 

~dV ~ A_ (l + £) ' A~ 

and A and i? vanish at infinity. 

Note that t, and r* are not well behaved functions at the horizon r = r + . 
As a consequence, the coordinate systems (t,r,9,<p) and (t,r*,9,<fi) only cover 

int (7£). Observe also that the two coordinate systems (t, r, 0, </>^ and (i*, r, 6*, </>) 
are identical for r > r cut . 

2.1.4 The Kerr-AdS metric for fixed (a,M,l) 

We may now introduce the Kerr-AdS metric as the unique smooth extension to 
1Z of the tensor given in the Boycr-Lindquist chart by: 



E , S , A e (r 2 +a 2 ) -A_a 2 sin 2 _ 
gKAds = ^-dr 2 + —dO 2 + i sin' 



A e (r 2 + a 2 ) - A_ . 2/1 , 7 , A_ - A e a 2 sin 2 9 , 9 
-2— ^A asin 2 6» d^dt dt 2 (9) 



where A_ is defined by Q and 



,2 



r 2 + a 2 cos 2 6», A e = 1 - — cos 2 9, E = 1 - ^ . (10) 



a 

p ~™ - - . p 

See [55] for explicit expressions for the inverse of ©. That the tensor © 
indeed extends to a smooth metric on 1Z is clear from expressing the metric 
in (t* ,r,9,4>) coordinates, which is carried out explicitly in [55]. Note that for 
a = 0, the metric (j9]) reduces to the well-known Schwarzschild-AdS spacetime 

/ 2M r 2 \ , 2 / 2M r 2 \ _1 , , 2 , 
ff = -^l-_ + _jdt 2 +^l-_ + _j dr 2 +r 2 ^ s2 , 

where dcrg2 is the standard metric on the unit sphere. 
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2.2 The norms 

Let d and V denote the induced metric and the covariant derivative on the 
spheres S 2 * r of constant t* and r in 1Z. 

We write \^...^ijj\ 2 = $ AA ' ...<f BB ' ^ A ...^ B ^ A , B ,ij) to denote the in- 
duced norms on these spheres. 

We denote by fij with i = 1,2,3 the standard basis of angular momentum 
operators on the unit sphere in 6, <p coordinates. 

With these conventions, we define the energy densities 



ei M = ^ (d t ^f + r 2 (d r i'f + |^| 2 + ^ , 

e 2 [i>] = ei [</>] + ei [d t *V>] + r 4 (d r c^) 2 + r 2 |d r ^f 



e2 [V'] = e2 [V 7 ] + 



f'i 



(11) 



Similarly, we define the following energy norms for the scalar field ijj, cf. [25] : 







r s ip 2 r 2 dr sin #d 


9d</> 










> (<w) 2 + 


iy^l 2 + V 2 


r 2 dr sin 




= u 


2 







r 4 (d r d r 4>) 2 + r 2 \fd r *P\ 2 + \ffi>\< 



r dr sin 9d9d(f> . 



Note in particular the relations 



/S t * 



eiM r 2 drsin^0= jlVllL.o _ 1 + H^^lllo, 



(12) 



P2 



' S t * 



r 2 dr sin OdOdcf) = \\Q,iip\\ 2 „ i,o rv > + |<9t*?A|| i-i-o , v * 

a AdS\ t*> AdS^ <•* > 

+\\nv° ( St *) + ii^ 9 '*vCo,- 2(Et (is) 



Higher order norms may be defined similarly. We denote by H AdS (Ylt*), the 
space of functions V 1 such that V*-0 S i 2 oc (S t *) for i = 0, fc and ||-0|| 2 fc s ,„ < 

oo. We denote by CH Ads the set of functions ijj defined on TZ, such that 
D C"(R t ,;H k A - q s ' S H^)) 

q=Q,..,k 

where Sfc = —2, Sfc_i = and Sj = s for j = 0, k — 2. 

When s = 0, we will feel free to drop the s in the notation, i.e. H Ads := H AdS 
and CH\ ds := CH^ dS . 
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2.3 A Final Remark 



In [26 , the coordinates t and <j> in © are defined with the x( r ) °f © being 
globally equal to 1. Here, for convenience in the subsequent analysis (which 
happens mostly away from the horizon, in r > r max ), we have altered these 
coordinates away from the horizon to agree with the Boyer-Lindquist coordi- 
nates. Note that these two coordinate systems are equivalent in the sense that 
the statement "||V'I! s ,^ decays logarithmically in t*" is independent of 

whether the coordinates (and E t * -slices) of [55] or the cut-off coordinates © 
are used. 



3 Separation of variables and reduced equations 
3.1 The (modified) oblate spheroidal harmonics 

For each £ € K, define the unbounded L 2 (sin 6 dOdifrj -self- adjoint operator 

P§2 (£) with domain being H 2 (S 2 ) -complex valued functions (see Section 7 of 
[15] for a more detailed discussion) as 

-PviQf = ^%(A«sin^/) + ^-^2- Q d V 
sin A e sin 9 f 



»|-cos* 0/-2z£— ^< 



We also define the operator P a , which is equal to 

flp,«(0 = { p (0 + ^ a 2 cos2 , ifo < . ( 15 ) 

For Z — > oo the operator Pga(£) reduces to an oblate spheroidal operator on § 2 
as considered for instance in [TS] . If also £ = a = we retrieve the Laplacian on 
the round sphere. Both operators Pg2^) and P§2 a (£) have discrete spectrum. 
The eigenvalues and eigenvector of P§2 Q (£) will be denoted as follows: 

P§ 2 a (£) bas eigenvalues A m £ (£) with eigenfunctions S m £ (£, cos0) e 1 " 1 ^. 

Later we will restrict attention to axisymmetric solutions and hence to the 
eigenvalues Xoe of the operators 

-Pe. a {e) :=-ifc,a(£) • (16) 

m— 

In fact, it will be convenient (in view of their manifest positivity) to work with 
the eigenvalues of the shifted operators (for a < 0) 

-Pe. a (e 2 ) - C 2 =— Q de (A e sm9d e (■))- ^a 2 cos 2 9 - ^a 2 £ 2 (17) 
v ' sine/ l z Aq 

and with the cos 2 replaced by a sin 2 in the second term in case that a < 0. 
In view of these considerations, we shall denote the eigenvalues of the operator 
Pe,a (£ 2 ) + £ 2 by Hi (£ 2 ). Note that by min-max and comparison with the 
spherical Laplacian (see Lemma 5.1 of |26j ) in the axisymmetric case we have 

lit (£ 2 ) > M£ (0) >Z£ (£+!)> c a ji (I + 1) . (18) 
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3.2 The separation of variables 

In this section, we present the reduced equations obtained after separation of 
variables. For the construction of quasimodcs, it would be sufficient to start 
directly from the reduced equations. However, to explain their relation with 
([T]), we will instead derive them from the Klein- Gordon equation ((T|). Thus, in 
this section ip will denote any regular solution to ([1]). 
Let us introduce the time-Fourier transforirQ 

1ft (t, r , = -L f e-^i, (u, r, 9, du . (19) 



We decompose the i\i of equation f|19[) as 

(aw) 
ml 



V (u, r, 6, </>] = (r) 5 m £ (aw, cos 9) e" n t> , (20) 

where S m g are the modified spheroidal harmonics introduced precedently and 
(V (r) = -=/ dt / dcre*"" S mt (auj,cos0)e- im4 ' ip(t,r,6,<f> , 



with d<7 = sin 6d6d<j>. 

After the rcnormalization 



S 2 (t,r) 



(aw) 



tfi? ) (r) = (r a + o a )»(^) (r), (21) 
we finally obtain from the equation 

(r)]"+(^-V^(r))u^=0, (22) 



where the potential (r) is defined as 



^ W - VIZ] (r) + (r) + V a (r) , (23) 



where 

r (au) 



o r 2 5^ + 3r 2 (l + «r -4Mr + a 2 

v +,ml\ r ) / 2 , 2^ C 2 , 2^3 

(H + a z ) (r z + a- 2 ) 

= ( r 2 +a 2 )~* (vWa 2 )", (24) 
A_ (A m£ + uj 2 a 2 ) - E 2 a 2 m 2 - 2mwaS (A_ - (r 2 + a 2 )) 



v 0,mt V ) 



(r 2 + a 2 ) 2 



V a (r) = - ~ A " 2 (r 2 + 9( a )a 2 ) , (25) 
1 (H + ar) 

with 8 (a) = 1 for a > and 9 (a) = if a < (recall that the \ m t also 
depend on a through l|15p). Note that V+ grows like =jp- near infinity while the 
Vo-part remains bounded. 



7 Note that in general iji is not an L 2 function in time and therefore. i/> is defined only as 
a tempered distribution. Since here we arc merely trying to justify the origin of the reduced 
equations, it will be sufficient to understand all computations formally. 
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3.3 The axisymmetric reduced equations 

We now look at the axisymmetric case, i.e. we consider the above equations 
under the assumption that ip is independent of the azimuthal variable <fi. The 
reduced equations are then obtained by setting m — in the decomposed equa- 
tions. Hence, we will consider the following set of equations: 

u[T ] (r)] " + (- 2 - W) (0 = 0, (26) 

where the potential (r) is defined as 

V { r } (0 = V jmik (r) + V mass (r) + V a (r) • W (cjV) , (27) 

where 

3r 2 5^ + 3r 2 f 1 + - 4Mr + a 2 

^•«nfc(r) = -A 2 — -r + A_ 



(r 2 + a 2 ) (r 2 + a 2 ) 

A-r 2 a A_ 2 
9 (a) a 2 , 



Vmass 



2 - a A_r 2 

Z 2 ( r 2 + a 2)2 ' 



K(r) = r 2 +~ 2 ^ (28) 

(r 2 + a 2 ) 

Here we rearranged the terms in the different potentials so that V mass = 
corresponds to the conformal case a = 2. In particular, we have \Vj un k (r) | < 
Cm i a , 2^~2i5 ^; ^ct an d hence that K-« n fc is uniformly bounded. 
The properties of are summarized in the following lemma. 

Lemma 3.1 (Properties of V a ). For all \a\ < I, the potential V a enjoys the 
following properties: 

• V c {r*) -> ^ as r* -> tt/2. 

• VJj /ias a unique local and global maximum V max at m [ri,7r/2). 
^Uso, V a is monotonically decreasing in [r* nax ,Tt/2). 

• V max = V a {r* max ) > I + 



Remark 3.1. In particular, for any < oq < I, and for all \a\ < oq, the size of 
the interval [1/7 2 , Kr(^ a x)] * s uniformly (in a) bounded from below by a strictly 
positive uniform constant. 

Proof. The first claim can be trivially checked. For the second and third claims, 
let us write 

A (r 2 + a 2 ) (£ + £ + z) -2Mr 
V a (r) = j = —^-L 

(r 2 +a 2 ) 2 (r 2 + a 2 ) 2 

1 5 (r 2 + a 2 ) - 2Mr 
P + (r 2 +a 2 ) 2 ' 



(29) 
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The r derivative of V„ is 
d T Vo (r) = ^~ 



2 (2Sr - 2M) - 4r (r 2 + a 2 ) (S (r 2 + a 2 ) - 2Mr) 



(r 2 + a 2 ) 4 
(2Hr - 2M) - 4r (H (r 2 + a 2 ) - 2Mr) 

r2 J- ^2^,3 



(r 2 + a 2 

-2Sr 3 + 6Mr 2 - 2Sa 2 r - 2M a 2 

r2 J- „2^3 



(30) 

(r z + a z 

If a = 0, then one has as usual that the only zero in [r+,oo) is at r = 3M. 
Assume a ^= 0. Observe that the derivative is positive near negative infinity, 
negative at r = 0, positive at r = r + and negative near infinity. This tells us 
that there are three real roots for d r V a . The one of interest to us is the one 
corresponding to the (unique) maximum of V a in the interval [r+,oo). Define 
r guess = tt". We claim r + < r guess < oo. This follows from r + < 2M, which is 

in turn a consequence of the fact that A_ = r (r — 2M) 

positive for r > 2M. Finally, at r guess = we have from ([29]) 

, i m + ^a 2 1 

V a (r guesa ) - ^ + +a 2) 2> p- 

V guess ^ " I 

□ 



4 Bound states 

As proven in our [26] . there exist no periodic solutions of the massive wave 
equation on Kerr-AdS. In this section, we will introduce an additional boundary, 
located at r = r max (the location of the top of the potential V a , as defined 
in Lemma 13. ip , enabling us to construct periodic solutions whose associated 
energies lie below the top of V a . 

In order to avoid confusion between the mode number £ and the real num- 
ber I determining the cosmological constant, we introduce the semi-classical 
parameter h > by defining 

h~ 2 (e,uj) := \ m {a 2 uj 2 ) + a 2 uj 2 = fi e (a 2 w 2 ) , (31) 

as well as the shorthand 

h Q 2 = h- 2 (£, 0) = fi e (0) . (32) 

In the rest of this paper, both notations, ui or Uh, will be used when we want 
to make explicit that a solution u depends on h or I. 

Having in mind a semi-classical type analysis with semi-classical parameter 
h, we then rewrite (|2"6l) as 

Non-linear Eigenvalue Problem: 

h 2 u 2 u = P (h) u := -h 2 u" + V a u + h 2 (V Junk + V mass )u , 

/■tt/2 (33) 

with boundary cond. u (r* nax ) = and / [\u'\ 2 + r 2 u 2 ] dr* < oo . 
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Remark 4.1. We will refer to the above boundary conditions as Dirichlet con- 
ditions. Note that they imply that u (tt/2) = 0. See Remark\4-2\ below. 



Unless a = 0, (|33p is a non-linear eigenvalue problem. Indeed, a solution to 
the eigenvalue problem 

ku = P (h) u (34) 

with Dirichlet boundary conditions is a solution of (f3"3"]l if and only if n = h 2 Lu 2 . 
In case a = 0, (|3"3"|) reduces to the linear problem ([M]) since the P (h) operator 
becomes independent of u), and, thefore, given any solution to (|34l) . one simply 
obtains a solution to (|3"3"f by defining ui 2 = h~ 2 n. 

What we would like to prove is that - given fixed parameters M, a, I, a 
- we can find, for any sufficiently small h (or cquivalently sufficiently large £), 
an uj 2 such that (|33[) is solved for some U£, and to control the size of this w| 
(cf. Proposition 14.21 below - ) . 

In order to understand (|33D . it will be useful to first look at the following 

Linear Eigenvalue Problem: 

h\io 2 u = P b a S e (ho) u := -h\u" + V a u + hl(V 3unk + V mass )u , 
with boundary cond. u (r* nax ) = and / [\u'\ 2 + r 2 u 2 ] dr* < oo . 

J r^ nax 

As explained above, ()35[) can be seen as a linear eigenvalue problem because ho 
and therefore Phase (ho) depends only on I (but not on oj). 

Remark 4.2. By Proposition 4 of J30\j . \35)) is a well-posed eigenvalue problem. 
This is a non-trivial statement because the potential V maS s is unbounded on the 
domain (r* nax ,n/2) unless a = 2. The condition f™/ 2 [\u'\ 2 + r 2 u 2 ] dr* < oo 

implies that ip = \/ r 2 + a 2 u ■ Sot (0) € H\ ds (in particular r 1 / 2 ~ e u (tt/2) = 
for any e > 0) and ensures the existence of a positive discrete spectrum with 
eigen] 'unctions in the energy spaced 



To prove that eigenvalues h^uj for (|35|) exist in a suitable range, we will 
perform a semi-classical type analysis for a semi-classical operator whose prin- 
cipal part should be —ti^u" + V a u. Since h^Vjunk is controlled by h^V^, this 
term will be lower order and hence negligible. On the other hand, unless one 
considers the conformal case a = 2, the a-priori lower order (in powers of ho) 
potential term h 2 ■ V ma ss is unbounded near r* = tt/2, so that some care (a 
Hardy inequality) is required. 

Observe finally that if we set a = in V a , Vj un k and V mass in (|35p , then 

Pbale (h ) U = K-U (36) 

would be precisely the eigenvalue problem one needs to study for Schwarzschild- 
AdS. In any case, in the next section we will establish the existence of eigenvalues 
k = h\u) 2 of the more general eigenvalue problem (|3"5"j) as the latter is easier to 
connect to the full problem (|33D . 



8 Using the twisted derivatives of 1301 . i.e. writing u" + VmassU = r n (r 2n (r n u)') ' + Vtwist ' 
u with n = 1/2 — l/2\/9 — 4o and hence Vtwist uniformly bounded, one could generalize the 
construction of the paper to other boundary conditions. 
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4.1 Weyl's law for the linear eigenvalue problem ( 135 h 

For the purpose of this section, it will be convenient to introduce the following 
notation. For all c < d, we define P^ff (c, d) as the following eigenvalue problem 

Pbase(ho) U= KU, 

with Dirichlet boundary conditions w(c) = u(d) = and - in case that d = tt/2 
- the condition dr* [\u'\ 2 + r 2 ■ u 2 ] < oo (cf. Remark IP)) , where Pbase{h ) 
is the operator defined by (f3"5| . Similarly, we will write P^ e (c,d) for the 
Neumann problem (which will never be considered with d = 7r/2). Finally, 
we write P^fff (c, d) for Neumann boundary at c and Dirichlet at d. In the 
latter case we again impose J dr* [\u'\ 2 + r 2 ■ u 2 ] < oo in case that d = tt/2 
to ensure that all eigenfunctions live in the energy space, cf. Remark l4.2l Note 
that P|) a i 5 e (r^ lax , 7r/2) is precisely the linear eigenvalue problem (|35|) . 
The aim of this section is to establish the following proposition. 

Proposition 4.1. Let a < f , M > and \a\ < I be fixed, and E € (4, V max ) be 
given. Then, for any 8 > such that [E — 5, E + S] C (jjr, V m ax), there exists an 
Hq > ; such that the following statement holds. For any < ho < Hq, there 
exists a smooth solution Uf lQ of the eigenvalue problem P^ e (r^ ax , 7r/2), with 
corresponding eigenvalue n lying in [E — 6, E + 6] . In particular, there exists a 
sequence ((/io) n , W(/ l0 ) ll ) such that the associated eigenvalues K((ho) n ) E as 
(h )„ -> 0. 

In the rest of this section, (a, M, I, a) are fixed parameters satisfying the 
assumptions of the proposition. 

We shall in fact prove in this section a stronger result than Proposition ^. 1[ 
namely a version of Weyl's law adapted to our problem. This is the statement 
of Lemma 14.21 from which Proposition 14.11 immediately follows. The proof of 
Lemma 14.21 in turn requires the following auxiliary lemma, which ensures non- 
existence of eigenvalues below a certain threshold. 

Lemma 4.1. Let E > be given. Then there exists an Hq > so that for all 
< ho < Ho, there exists a r* K (E, ho, a) such that the problems P 1 ^^ (r* K , tt / 2) 
and P^f ) e (r* < ,ir/2) have no solutions with k < E. Moreover, 

fmax <r K < j- ■ E, (37) 
ho 

where C depends only M,a,l and a. 

Proof. Assume there was a solution u of Pj)i) e {r* K , n/2) or P^i) e {r* K , tt/2) with 
k < E. Then we would have 

tt/2 

dr* [h 2 \u'\ 2 + (V a + hlV junk + h 2 V mass - E) \u\ 2 ] < (38) 

K 

for this u. On the other hand, the Hardy inequality 

tt/2 a ,tt/2 

dr* \u\ 2 < 4Z 2 / dr*\u'\ 2 (39) 
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proven in our |26j holds for u. This implies that 

< 0. 



tt/2 

dr* 



l o^^J~^2 + v <? + h>lVj Un k + hlV mass - E ) \u\ 2 



The dominant term in the integrand near infinity is fig (| — a) ji which is pos- 
itive, while all other terms remain bounded. Hence by choosing r* K sufficiently 
large {tk > 7— • E for some constant C) we obtain a contradiction as the round 
bracket in the integrand eventually becomes positive. □ 

Consider now the eigenvalue problem P D a ff (r* lax ,ir/2) and fix an energy 
level £ £ (jj,V max ) . Lemma l4~Tl produces an r* K (£,h ,a), to which we associate 
the phase-space volume 

V £ , ho , a = Vol{ (£, r*) £ R x [r* max y K ] \ ? + V a + h 2 V mass + h 2 V jank < £} 



dr*y/£-V(r*)- X{ v(r*)<s}- (40) 

Note that for fixed £, this expression converges uniformly in ho as ho — > 0. This 
is already immediate for a < 2: V (r*) is then bounded below and hence the 
integrand itself is obviously uniformly bounded in ho- For 2 < a < 9/4, the 
integral (|4T)|) also converges uniformly in ho since 



: i 

dr*h a r <C / drh a - < Cho log 



/l 



goes to zero as /i — > 0. Here we have used the estimate ([57)1 on . 

Finally, to state and prove Weyl's law, we also introduce an expression for the 
phase space volume between two energy levels, say [E — S, E + S] C (j?,V max ): 

Q.E,a = lim VE+8.h Q ,a — Hm VE-5,h a ,a 
ho— >0 ho— >0 

= Vol{(£,r*) \ E-5<£ 2 + V a <E + S}. 

By an elementary computation, we have a lower bound Qe,u > Ce.m.i,u • 5 for 
a constant independent of ho- 

Lemma 4.2. Consider the eigenvalue problem P^f> e (rj,^, tt/2). Fw; an energy 
level V max > E > -p and prescribe a S > small. Then the number of eigen- 
values of Pj^f (C.,7r/2) Zj/m# m tfte interval [E — S,E + 5] C 
denoted N [E — 5, E + 5] , satisfies Weyl 's law 

N[E-6,E + 6\ ~ -L-Q Eia . (41) 
ho-yo Zirho 

Proof. Choose r* K [E + S, ho) such that by Lemma \A. II there are no eigenvalues 
below E + 5 of P h D a D e {r* K , tt/2) and P^ e {r* K , tt/2). We equipartition the domain 

[ r max7 r if 0^! ^o)] into intervals of length /3 = -^— k — SL - We then consider the 
following two comparison problems: 

• The Dirichlet problem Pjjff (r* K , tt/2) in conjunction with k Dirichlet 
problems indexed (i = 1, ...k) and arising in the following way: They 
are the problems P D a ff (r* max , r* max +/?), ..., P D a £ e (r* K - f3,r* K ) but with 
the potential replaced by a constant, which equals the maximum of the 
potential on the interval. 
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• The mixed problem P^rf ( r /o in conjunction with k Neumann prob- 
lems indexed P^~ (i = 1, ...k) and arising in the following way: They 
are the problems P b N % e (r* max , r* max + /?),..., P b N % e (r* K - p, r* K ) but with 
the potential replaced by a constant, which equals the minimum of the 
potential on the interval. 

We can estimate the number of eigenvalues of P^ff (r* nax , tt/2) below a thresh- 
old £ by 

k 

£ N< £ (P*+) + N< £ (P b D ^ (r* K> tt/2)) < 

i=l 

k 

N< £ (P b D T (r* max , tt/2)) < ^ iV< £ (i*") + 7V< £ (P^ e (r£, tt/2)) . 

By our choice of r* K , we have N< £ (i*g e (r* K , tt/2)) = N< £ (P#£ e (r* K , tt/2)) = 
for £ = P + <5. On the other hand, for each Pp~ and each P^j 7 , the number of 
eigenvalues can be estimated directly (as each problem can be solved explicitly). 
We have 



l+ ) 



; = 1 




The estimate for P l N is similar with the potential being replaced by VI and the 



number of eigenvalue in each cell being m ax (0, £ — V!) y £ — VI +1. 

To conclude, let us choose the number of cells k such that k(ho) tends to oo 
as ho goes to and moreover k(h) = o{l/ho). The sums converge as a Ricmann 
sum and the errors are then of order o(l/ho). Therefore we get 



g N * ^ ho Zo 2k [i dr * V^vW) ■ xnm<_ £ - (42) 

The statement of the Lemma then follows from 

N[E - 5,E + 5] = N< E+S - N< E -s 
using the previous formula with £ = E ± 6. □ 

4.2 Kerr-AdS 

In the last section we showed that for any fixed given parameters M > 0, \a\ < I, 
a. < |, the eigenvalue problem (|35p . Phase (ho) u = k-u with Dirichlet conditions, 
admits (lots of) eigenvalues k in the range E — 5 < k = h\io 2 < E + S, provided 
ho is chosen sufficiently small (i.e. £ large). 

As an immediate corollary, we obtain the existence of eigenvalues in the 
desired range for Schwarzschild-AdS, simply by setting a = 0, cf. ([Ml) . For the 
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Kerr-AdS case, we still need to relate the above result to the full problem, which 
we recall is the non-linear eigenvalue problem Q33[) given by 

P(h)u = n-u with K~Lo 2 h 2 

and boundary conditions u {r* nax ) — as well as f*/ 2 dr* \\u'\ 2 + r 2 u 2 ] < oo. 

To achieve this, consider for fixed \a\ < I, M > 0, a < |, the two-parameter 
family of linear eigenvalue problems 

Qe (b 2 ,uj 2 )u = A(b 2 ,uj 2 )u (43) 

for the operator 

Qe {b 2 ,uj 2 ) u := -u" + (V a ft e (b 2 a 2 uj 2 ) u + (V junk + V mass - a; 2 )) u , (44) 

complemented by the above boundary conditions. Here b 2 G [0, 1] is a dimen- 
sionless parameter and ui 2 G M + . Our goal is to show that for b = 1 there 
exists an u> 2 such that the above problem has a zero eigenvalue and to moreover 
suitably control the size of this w 2 . 

By the results of the previous section, we know that for 6 = there ex- 

ists, for any sufficiently large t, an us 2 e (satisfying E — S < ^jjy < E + S) 

such that Qe ^0,oJq£^ admits a zero eigenvalue. Moreover, this eigenvalue is 
non-degenerate by standard Sturm-Liouvillc theory Listing the eigenvalues of 
Qe ^0,Wq in ascending order, let us say that it is the (ne) th eigenvalue, which 
is zero. 

The strategy, now, is the following: We will show by an application of the 

2 

bJ 



implicit function theorem that for any b € [0, 1] we can find an oj 2 t such that 
the (ne) th eigenvalue of the operator Q (b 2 ,ui 2 A is zero. As a second step, we 



u b,e 

will provide a global estimate on the quotient — -, - b ' e „ r . For this last step an 

w( b a u b,t) 

important monotonicity will be exploited. 

Lemma 4.3. Fix parameters \a\ < I, M > and a < |. Suppose we are 
given parameters bo £ [0, 1] and uj 2 () e G IR + such that the eigenvalue of 

Ql [b\ , uj 2 o is zero. Then, there exists an e > such that 

1. for any b 2 G (max(0, 6§ — e), 6q + e) one can find an associated ui 2 ^ G IR + 
such that the (ne) th eigenvalue of Qe (b 2 ,uj 2 ^j is zero, 

2. ui 2 e changes differ entiably in b 2 G (max(0, 6q — e), 6q + e) and we have the 
estimate 

< < Cut* 

for some constant C > which is independent of bo (but may depend on 
I). 

3. The e > can be taken to be independent of bo (but may depend on I). 
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4- ui 2 e satisfies the estimate 



c- L < 



J b.t 



1) 



< c, 



for some c > depending only on the parameters a,l,M,a. 

Proof. The n th eigenvalue of Qi (b 2 , ui 2 , denoted A n (fe 2 , w 2 ) , moves smoothly 
in the parameters b 2 and w 2 and we have the formula 

A„ (b 2 ^ 2 ) = dr*V« <3f (& 2 V) V„ (&W) (45) 



for the eigenvalue, provided we normalize the associated eigenfunctions ip n (b 2 , uj 2 
by J*/ 2 dr*\ip n (6 2 ,w 2 ) | 2 = 1. By assumption, A„ (b 2 , 0J 2 o ,e) = 0- We compute 



duj 2 



Ji (^U - jH~ ■ (V ®><J ) ■ (47) 

The angular eigenvalue m (b 2 a 2 oj 2 ) is itself a smooth function of the two pa- 
rameters b 2 and ui 2 . We have the formula 

fit (b 2 a 2 Lu 2 ) = f dOsmO fa (b 2 ,oj 2 ) [P e , a (b 2 a 2 uo 2 ) + b 2 a 2 uj 2 ] fa (b 2 ,uj 2 ) (48) 



K<J = C dr^l (b 2 ,u 2 oA ) ■ (V a ■ (6§,<<) - l) , (46) 



tt/2 



9/^ 



for the eigenvalues provided we normalize the associated eigenfunctions fa (b 2 , uj 2 ) 
by f£d0 sin 6 \fa (b 2 ,uj 2 ) \ 2 . Recalling from (JI7J) that 



1 



P g , Q (b 2 aW) f + b 2 a 2 cu 2 f 

a / a . «o sm2 ^,2 2 2 , M 2 f COS 2 9 f 

-d g (Aesmedef) + — r — & 2 aV/ + ^a 2 ' 



(49) 



smff " • A e '' ' Z 2 

to be read with the upper (lower) line in case that a < (a > 0), we obtain 

^K<d=b 2 l 2 -T x and ^ (& 2 , , (50) 

where 



T T = / d0 sin 



sin 2 a 2 



2 , ,2 \ |2 



< sup 



sin 2 t 


?a 2 ~ 


/ 


[ A e 


I 2 _ 





£d6 S m6 \fa{b 2 ,u 2 bo>l ) | 2 <^. 



(51) 



where we have used the estimate supggjo.Trl S A e — 1; which can be easily 
checked using (|TU)) . Going back to (|4l)| and (|4T|) , the implicit function theorem 
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allows us to solve for uj 2 as a (smooth) function of b 2 locally near U)q,U! 2 o A, 
provided that the right-hand side of (|46[) is non-zero. To achieve this, note that 



1 - > 1 - bla 2 V n > 1 - a 2 V a 



du 2 



using the estimate obtained for from ((50|) and (|5T|l . On the other hand, 
using ([29) 

2 i — i 9 9 i — . 9 

2 a z cr „ 

a ^ = 77 + 7^2 =72+-- 77^772 <1 - CA ^- 

where CM,i,a is constant depending only the parameters M,l,a. It then follows 
that the right-hand side of (|46[) is bounded away from zero with the lower bound 
being independent of bo. This concludes the proof of the first item of the lemma. 

The implicit function theorem also provides a formula for the derivative of 
the function w 2 (& 2 ) just obtained, namely 



In view of (|50p. we immediately obtain an estimate of the form 

0<^(6o)<C MM <,, (53) 

establishing item [2] of the lemma. As a consequence, we also obtain a uniform 
bound on ui 2 e for all b € [0,1]. 

Observe on the other hand that the identity A„ (w 2 ,6 2 ) = implies (us- 
ing the Hardy inequality (|39[) as well as V a > 77) that w 2 £ > ^p/ie (b 2 uj 2 ^ > 
2fi^e (0) > CM,i,a£ + 1) for any 6 2 € [0, 1]. Hence the quantity w 2 £ will stay 
strictly away from zero and 

,,2 , ,2 

Cma ° - 7(7TT) - ^'-IcF+i) - CM,i '° ' (54) 

the first inequality following from the analysis of the linear eigenvalue problem 

(|3"5j) in Section 231 anc » the second from the uniform estimate (|53[) on .A* . 

Finally, the size of the e promised by the implicit function theorem is uniform 
in b 2 by the fact that A n (b 2 ,oj 2 ) is a smooth function on the compact set 



[0,1] x 



0,w 



2 

max, J 



, where ^ denotes an upper bound for uj£ 2 r This, together 

with the uniform estimate (|51[) allows uniform control of the errors arising in 
the implicit function theorem. □ 

In view of the fact that e is uniform in 6q , we can apply the implicit function 
theorem all the way from 60 = to b = 1. Note that e can depend on I. 
However, for each fixed £ it only takes finitely many applications of the implicit 
function theorem to reach 6=1. 
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To gain quantitative global control beyond ([54"]) on the behavior of co 2 e {b 2 ), 
let us look at the quotient 



E n (b 2 ) 



J b.( 



[9) 



with b 2 G [0, 1] 



(55) 



W (Ww 2 , (6 2 ) 

which by construction is the (ni) th eigenvalue of the semi-classical operator 
Qb (it) = E n (& 2 ) u with 

Qb (u) = -u" [fie (b 2 a 2 u> 2 tl )] + V a + [fit (6 2 a 2 w^)] 1 (V junk + V mass ) u . 

Recall that E is an energy level such that E £ (jy, V max ) and that, with the 
notation just introduced, E n (0) £ [E — S, E + 8) for some small 8 > 0. 

Lemma 4.4. For all 8' > 0, there exists a L, such that for all i > L, we have 

1 



r 2 



< E n (l) <E + S + 8'. 



(56) 



Proof. We first establish the upper bound. Using the Hardy inequality ([55)1 
together with (fT5|) (which implies that a := j^qj — ,2 \ > f° r am/ 



1 > b > 0), we can estimate 



E (b 2 a 2 



«r/2 








/ dr* 


[.( 


Qo- 


Qbj u 


T max 









> a 



> o 

tt/2 



tt/2 



dr*\ 



dr*\u\ 2 
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1 A_ 

4/ 2 r 2 + a 

/ 2 r 2 + a 2 



C M ,l,a\ U Y 



(57) 



since Vj un k is bounded uniformly. In view of \cr\ < ^^pjy, we conclude that 



Qo> Qb — 



holds for any 1 > b > 0. Hence, by min-max, we infer in particular that 
(independently of the parameters a, M, I and a) 



E„ (1) < E n (0) + 8' < E + 8 + 8' , 



(58) 



where 8' can be chosen arbitrarily small by choosing £ sufficiently large. 

For the lower bound, we will establish that J dr*u (Qb (u) uj > -p holds 

for any u £ H\ dS and I large. To prove the latter, note first that the Hardy 
inequality (|39p reduces the problem to showing that 



Va 



1 A. 



4 r 1 + a' 



~t~ Vjunk ~f~ Vmass 



> 



r- 



(59) 



The square bracket is manifestly positive in a region [r£,7r/2) for some large r* L 
close to 7r/2 (depending only on the parameters M, \a\ < I and a < §). Wc fix 
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this r* L and, in view of V a > A, have established (f5§)) in [r* Ll ir/2). In [r^ ax , r£], 
we have the global estimate 

1 Mr 

The second term on the right will dominate the term h 2 ■ (Vj un k) pointwise in 
[ r moi! r i] provided h is chosen small depending only on M, I, a, a. □ 

We summarize our results in the following proposition, which can be under- 
stood as the analogue of Proposition ^. II 

Proposition 4.2. Let a < %, M > 0, \a\ < I be fixed and E € (jT,Vm ax ) 
be given. Then, there exists an L > such that for any i > L, the following 
statement holds. There exists an ui 2 £ M + and a smooth solution ut of the 
axisymmetric reduced equation 

- u'l + \V a ■ fit (a 2 0jj) + V junk + V m ass] Ui = U}j • U t 

/•t/2 (61) 

satisfying ue (r^aa:) = and / [\u' e \ 2 + \ut\ 2 r 2 ~\ dr* < oo . 



Moreover, the uij satisfy the uniform estimates 

fit \a l ^i) 2 £ + lj 



5 Agmon estimates 

In this section, we recall the so-called Agmon estimates. These are (well known) 
exponential decay estimates for eigenfunctions for Schrodinger type operators, 
in the so-called forbidden regions. 

5.1 Energy inequalities 

The Agmon estimates will rely on the following identity 

Lemma 5.1 (Energy identity for conjugated operator). Let r* > r$. Let h > 

and let W , <f> be smooth real valued functions on [r*, r*J. For all smooth functions 
u defined on [r*,,r*], we have the identity 

W - 

\ \dr* J J ' J 
Wuj ue 2 idr* + J h~ 1 ^e 2 %2i% ^u-^ ] 

In particular, if u is real valued and vanishes at and r*, then, 

I e 2 i\u\ 2 | dr* = 







1 + h~ 2 Wu ) ue 2 %dr* . 



dr* 2 
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Moreover, the same idendity holds if W is not assumed smooth on [rg,r*] but 
only such that W\u\ 2 £ L l {r^,r\). By density, we may also replace smoothness 
of u and <f> by u G Hq [Vg , r\] and (f> is a Lipschitz function. 



Proof. This follows easily from the computations 



dr* 2 



and 



— + h W I (e7r u)ue~z dr* 



dr* 



-e h u 



h- 2 We 2 ^\u\ 2 dr* 



d . ± . d , $ 

— (ehu)ue" (rj + y^Ce* 
dr* J dr* 



u)ue*n (7"q) 



= r* 



-h 



— 1 d<j> du 



dr* dr 



-ue h 



- h 



dr dr 



dr- 



dr* 

2 



dr* 



J!] 



4 (&) l-l 2 - 2 

+ fj h- 1 ^-e 2 i2i%(uu r *)dr* 
-h- 1 (M\ v \^)(r t ) + h-^(M\u\^) 



□ 



5.2 The Agmon distance 

We will rely on the Agmon distance to establish our exponential decay esti- 
mates]! Given any energy level £ > 0, and a potential V = V(r*) (which may 
also depends on an parameter h), we define the Agmon distance d between r* 
and r\ as 



d = d(y-e) + (fi,rZ) 



X{V>S}(r*)(V(r*)-£) 1/2 dr* 



where X{v>£ } is the characteristic function of the set of r* satisfying V(r*) > £. 
In other words, d is the distance associated to the Agmon metric (V — £)+dr* 2 , 
where / + = min(0, /) for any function /. 

It is easily checked that d satisfies the triangular inequality and that 

\V r <d(r\r*)\ 2 <(V-£) + (r*). 



9 The Agmon distance is actually typically used to obtain optimal exponential decay esti- 
mates, see for instance 1201 . For the main purpose of this paper (the construction of quasi- 
modes), we could have used smooth cut-off constructions to prove slightly weaker exponential 
decay estimates. However, the Agmon distance (despite leading only to Lipschitz cut-offs) 
has a nice interpretation which is why we choose to use it here. 
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The distance to a set can also be denned as usual. In particular, we define 

d £ {r*):= mt d(r*,r* ), 

which measures the distance to the classical region. We have again 

\V r >d £ (r*)\ 2 <(V-£) + (r*). 
For a given small e € (0, 1) we define the two r*-regions 
17+ = fi+(£) :={r* | V{r*)>£ + e} 

and its complement 

Q- = Q~(S) := {r* | V(r*) <£ + (}■ 







^ c -t- e 

\ £ 















5.3 The main estimate 

We would like to apply Lemma 15.11 between r^ ax and n/2 for u a solution to 
the eigenvalue problem ([55)1 and for suitable 0. 



Lemma 5.2. Let u be a solution to the eigenvalue problem 
P(h)u for some n = h 2 u 2 . Define for any e € (0, 1) 

(f> K ,e ■= (1 - e)d K . 
Then, for all e sufficiently small, u satisfies 

fir/2 



i.e. k ■ u 



(63) 



d *m 
— — e h u 
dr* 



dx 



e- / e- 



\u\ 2 dr* 



< D(n + e)e 



2a(e)/hl 



(64) 



where a(e) = sup^- d K and D > is a constant depending only on the parame- 
ters a, M, I and a. 



Remark 5.1. Note 



a K (e) = supd K -> 



as £ -> 0, uniformly in h ( and K ) for h sufficiently small. In view of the 
exponential weight in the second term on the left, the estimate {64\j quantifies 
that u is exponentially small in the forbidden region, provided we can show a 
uniform lower bound for (f> K _ e in a suitable subset of Qf . This will be achieved 
in Lemma \5.3\ . 
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Proof. Applying Lemma between r^^ and 7r/2, we get 



TV /2 



dr 
K-V 



e h u 



dr* 



dr* 



at 



e 2ft lit! 



r/c? 



e 2 *|w| 2 dr* 



2 dr*. 



In view of our choice 4> = </>« e , we have in f2+ the estimate 



dr* 



> {l-{l-ef){V-K)> 



(65) 



(66) 



for e sufficiently small, which we will use to estimate the left hand side of (|65j). 

For the right-hand side of (|55|) . we note that if V > (which occurs if a < 2), 
then we immediately obtain 









2' 


/ 


k-V + 






IQ7 




dr* 





e 2 ~^\u\ 2 dr* < (n + e)e 



2a(e)/h\ 



lL 2 (r^ ax ,7r/2) J 



(67) 



so that combining (|6"6) and (j6"T|) yields 

To obtain also in the case a > 2 (for which we have V(r*) — > -co as 
r* — > 7r/2), we need once again to appeal to a Hardy- type inequality to absorb 
the error by the derivative term on the left-hand side of (j65|) . 



This we do as follows. Recall that V = V a + h 2 (Vj un k + V mass ), and that 



the unbounded term is h 2 V ma 

A_r 2 A_ 



h 22 ^ , < for a > 2. 



Note that 



V~ (r 2 +a 2 ) 2 
A a 2 

f „ 2 Z„2\i ■ The second term is bounded (and 



(r 2 +a 2 ) 2 (r 2 +a 2 ) (r 2 +a 2 ) 

in fact will contribute with the right sign if a > 2) so its contribution can be 
treated as before. Thus, we only need to estimate f n - ( r ^~ 2 ) e 2 ~^~ \u\ 2 dr* . 
By Lemma 7.2 of [55] (cf. we have for any function v in -ffp (r^ ax , 7r/2) 



OT /2 A _ 



(r 2 + a 2 ) 



\v\ 2 dr* < 4l 2 





dv 


Ir* 


dr* 



dr* for any R* > rj^ 



(68) 



Applying the above Hardy inequality to v = e^ K -'^ h u, we obtain that there exists 
a uniform constant C > such that 



tt/2 



h 2 (9/4-a) 



d 

— — e u 

dr* 



dr* 



V 



dr* 



e 2 -^-\u\ 2 dr* 



<C(n + e)e 2 ^ h \\u\\ 2 L2{r ^ /2) , 
i.e. there exists a constant Z? > (which degenerates as a — > 9/4) such that 



tt/2 



dr 



-e h u 



dr* 



dr* 



e 2 -^\u\ 2 dr* 



< D( K + e)e 2 ^^\\u\\ 2 L2{r ^ /2) . (69) 
This estimate, when combined with yields again (|64[) from (|6"5"|) . □ 
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5.4 Application of the main estimate 

Before we can exploit we need the following Lemma, which quantifies the 
size of the forbidden region for a given energy level. 

Lemma 5.3. Let E G (p-, Vmax) and suppose that k € (j?>E + <5] f or some 
S > smc/i that E + 5 < Vi nax . TTien £/iere exists a 6' > and a C > 0, 6<rf/i 
constants being independent of h, such that V a — k > 2C, in [r max ,r max + 5], 
for all k£[E- 5,E + 5}. 

Proof. This is a simple consequence of the continuity of V a at r* nax . □ 

In view of the full potential being V = V a + h 2 (Vj un k + V mass ) we also 
obtain 

Corollary 5.1. For h sufficiently small (depending only on M, I and a) we 

have V — k > C in [r max , r max + 5'] for all k G [E — 5, E + 5] with both 5' and 
C depending only on M , I and a. 

With E g (7*r,Vmax) given, we now fix 5' > and C > as promised by 



Lemma 15.31 This implies that 



> c M ,i.a in [r. 



max ; ' max 



5'] uniformly in e 



(the constant CM,i,a being of size C • 8'). Next we fix e > sufficiently small 
so that in particular a(e) < cm, 2, a/2. We finally conclude from that there 
exists a C > (independent of h) such that 



\u\ 2 dr* < Ce- c/h \\u\\ 



(70) 



Turning to the derivative term on the left of (p^)) . we also have 



/: 

J r* 





d(j) K , € 


(h 2 


dr* 



2 d4> K>IL du 
h dr* dr* 



du 



dr* 



3 2a(e)//i| 



li 2 (r 



dr* < 

^x^/2) 



The \u\ 2 term in the above integral can be ignored since it has the right sign. The 
crossterm can be absorbed using ([70)) and \ of the derivative-term. Therefore, 









fr 







dr* <C^ 2 e- 5 /' l || U || 2 r 



L 2 (r* lax ,-n-/2) 




Summarizing these decay estimates, we have proven: 
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Lemma 5.4. Let E 6 (jjj-, V max ) be fixed and let 8 be sufficiently small so that 
[E — S, E + 5} C (jz,V niax ). Then, there exists constants D,S' > depending 
only on the parameters M , I, a and a, such that the sequence of eigenfunctions 
[ u i\t>L arising from Proposition ^. 6 ^ satisfies the estimate 

+ \uA dr* < De- D / h \\u e \\ 2 LHr ^ w/2) , 

where h = (u^ (a 2 w 2 )) L/l and ui£ are as in Proposition \4-%\ 

We remark that by reusing once again the equation, we can obtain such an 
exponential decay estimates on all higher order derivatives, with the constants 
in the above lemma depending on the order of commutation. 




6 The construction of quasimodes 

By now we have established the existence (Proposition I4.2[) of a sequence of 
functions [ui] such that for each £ the corresponding ug solves 

io\h 2 ui = P(h)ue, 

— 1/2 

where h = (/i^(a 2 o;|))) — )• as I ~ > oo, and such that these ue obey the 
estimate of Lemma 15.41 with some constants D,5' > independent of h (or 
equivalcntly £). 

Let now x be a smooth function such that \ = 1 on [ r max + ; 7r /2] and 
X = on (— oo,r* lax ]. We then define ijje(t,r,0,(j)) as 

Mt, r, 6, $) = e^'x^MXr 2 + a 2 )" 1 ' Vr*(r))S«> (*)• ( 71 ) 

Remark 6.1. As defined above, the ipi are complex functions, but of course, 
we could have worked below with Re(ipi) or Im(ipi). 

The next Lemma shows that the i^i satisfy the Klein-Gordon equation up 
to an exponentially small error: 

Lemma 6.1. For each I and each k > 0, ipg £ CH\ dS . Moreover, there exists 
L > such that we have the following estimates. For all k > 0, there exists a 
Ck > such that for all I > L, for all t > to, 

Finally, all the H\ ds norms of ip£ and of its time derivatives on each S t are 
constant in t. 

Proof. Note that by standard elliptic estimates, any ui is smooth on (r* nax , tt/2). 
Thus, as far as the regularity of ip( is concerned, it is sufficient to check that ipi 
and its derivative decay sufficiently fast near r = oo, which is easy and therefore 
omitted. 
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Moreover, in view of our construction, we have O g ip£ = in [r + ,r max ] U 
[ r ( r max + $')■> °°]- Hence, the error is supported in a bounded strip in which we 
have the following naive estimate: For all (t, r, 9, </>) with r € [r max , r max + 5'], 

\a g ^ + < (^ 2 kl + Kl + Kl + /^ 2 M) aeo(fl), 

which gives the required estimate for k = after integration, using the Agmon 
estimates of the previous section and the equation satisfied by U£ in order to 
estimate u'[. For higher k, it suffices to commute the equation and to use the 
equation for U£ every time two radial derivatives occur, or the equation for Seo 
every time angular derivatives occur. □ 

Note that we finally proved Theorem 11.21 Indeed, the ipe arc of the form 
claimed in[T] by construction of (|TTj) . The estimate on the ue in[H was obtained 
as part of Proposition 14.21 The error-estimate |3l ist the statement of Lemma 
16.11 while the localization properties |U and O are obvious from (jTTj) itself. 



7 Proof of Corollary 11.11 

In this section, we prove Corollary 11.11 Given the quasimodes, the proof is 
standard, but we include it for the paper to be self contained. 

Let us therefore fix a Kerr-AdS spacetime such that the assumption of Corol- 
lary [O] are satisfied and also a Klein-Gordon mass a < 9/4. For convenience, 
we set tg = 0. Recall also that t* = t in r > r max . 

We shall consider solutions ip to homogeneous and inhomogencous Klein- 
Gordon equations with initial data i/'lst an< i dtip\^ t given on slices of constant 
t. We shall avoid completely issues rcgarging the facts that dt is not always 
timclikc and that the coordinate t breaks down at the horizon by considering 
only axisymmctric data which is compactly supported away from the horizon. 

Thus, given any (,s6l and given any smooth, axisymmctric initial data set 
w = Gp, Tip) , whose support is bounded away from the horizon and which decays 
sufficiently fast near infinity, we will denote by P(t, s)w the unique solution at 
time t of the homogenous problem 



P 

dip 
dt 



ijj = 0, 



1>, 
J\p. 



Given a smooth axisymmctric function F defined on 1Z, compactly supported 
in r away from the horizon and infinity, we can consider the inhomogeneous 
problem 



dip 



F 

J\p. 



For regular data as above, this problem is well-posed in CH\ dS and we shall 
denote its solution by ^p(i), suppressing the dependence on r and the angular 
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variables. If the data is axisymmetric, then ipp will be axisymmetric and writing 
v(s) = (0, F(s)(g tt )~ 1 ), ipF(t) is given by the Duhamel formula 

ip F (t) =P (t,0)w+ [ P{t lS )v{s)ds. 
Jo 

We now consider the family of ipi given by Theorem 11.21 For each I , ip£ 
provides an initial data set wg = (^pe(t = 0), ^f(t = 0)j for ([T]) on the slice 
t = 0. Moreover, ip£ satisfy the inhomogeneous Klein-Gordon equation 



for some F# satisfying \\FA\ u k,-2 < C^e HV^IIfr 

AdS AdS 

Let ip£ denote the solution of the homogeneous problem associated with the 
same initial data we, i.e. ip£ = P(t,0)w£. From Duhamcl's formula, we then get 



H , 



,(E t n{r>r max }) 



< 



t sup ||-P(*,s)(0,F^)(s)|| iJ i (is(Etn{r > rmax}) 

s€[0,t] 



< *C||fi|| lr o.- a(Eo) 



< tCe ||^||ifO ds (Son{r>r max })7 



tCe 



-ct 



; (E n{r>r max }) 



(72) 



where we have used the boundedness statement of Theorem 11.11 to bound 
||P(i, s)(0, ^)(s)||ffO dS (s t n{r>r max }) m terms of the data, as well as Lemma RTT1 

In particular, since the norms of rp£ are time invariant, for any t < 
reverse triangle inequality and f|T2[) yield 



2C 



the 



ei 



'E t n{r>r max } 



r dr sin 



> 



,J t\\H\ dg (E t n{r>r mx }) 



,(E n{r>r max }) 



> 



2( 



l fi i^llHi dg (Eon{r>r mM }) + I \ 9 t* Ipt I \ll\ is (S n{r>r max }) 

r 2 <ir sin OdOdcf) 



> 



2( 



(73) 



Here we have used - in the step from the second to the third line - that the 
data for ip£ is frequency localized, which allows to exchange angular and time 
derivatives with powers of £ using the second item of Theorem 11.21 and radial 
derivatives by angular and time derivatives using the wave equation the ip£ 
satisfy. From the third to the fourth line we exploited the fact that the data is 
localized in r > r max . The constant c depends only on the parameters M, I, a 
and a. 



Finally, setting t£ 



2C 



obtain from (|73|) a family (tt,ipi) such that, for 



large enough, Q ijj£ (tf) > C > holds for any £, which proves the Corollary. 
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A The Improved boundedness statement 



The boundedness statement at the if 2 -level proven in [35J 121] is the estimate ^ 
for the t2 [V'J-based energies, cf. Section |2~21 It is remarked in [25] that stronger 
norms can be shown to be uniformly bounded using commutation by angular 
momentum operators leading to the statement ([3]). Since the latter statement 
has been used in this paper and also in [26] . we provide here a sketch of the 
proof of this well-known (but absent from the literature) argument. We define 
the energies 



d [ip] (f) r 2 sin 9drd0d(/) 



(74) 



and with the obvious replacement, E2 [ip] (t*) and E2 [ip] (t*). Recall that uni- 
form boudedness for the E 2 [tp] energy is derived, in addition to known tech- 
niques near the horizon (cf. the red-shift vector field), by commuting the Klein- 
Gordon equation with dt (which yields ([2]) with tp replaced by dtip) followed by 
elliptic estimates on spacelike slices, which control the H\ dS norm. 

Let us sketch how to prove boundedness for the E2 [tp] energy. If we com- 
mute the Klein-Gordon equation with angular momentum operators we obtain 



n g (Q i i>) + p(Pii>) 



7T ■ VuV^V ■ 



2V C 



1 afi 



(75) 



with ^*'7T the (non-vanishing in Kerr!) deformation tensor of The right 
hand side decays suitably in r but not in t. More precisely, in view of the fact 
that there is no integrated decay estimate available, we cannot close the basic 
energy estimate on its own. Let us instead commute with localized angular 
momentum operators Qi = \ (r) f2^ where \ ( r ) is equal to 1 for r > 2R and 
equal to zero for r < R. Applying the energy estimate for the vectorficld dt* , 
we can derive 



1 1 ^ii>t 1 1 _f/°.- s 2 (E t , n{r>2R}) 



10,; 



\H\, S (£ t *n{r>2fl}) 



< CM,l,a,ct 

+ to - n) 



11^^11^0,-2^ ) + 11^11^(13,3) 



e sup E 2 [tp] (r) + e • E 2 [ip] (0) 

tE(ti,t 2 ) 



(76) 



where e can be made small by choosing R large. The last term arises from the 
spacetime error-term which decays strongly in r. 

The idea is to combine this with an integrated decay estimate for the Vlitp 
which loses linearly in r. Recall that if Oty + jpty = /, then we have the identity 



a 



fif t " t "Vi**V t **-25 t 6 V f .$Vi,$-/* 



(77) 



where a, b run over r, 9, <fi only. When integrating this identity (with ^ replaced 
by ttiip and / being the error arising from the commutation in (1751) ) with the 
usual spacetime volume we observe that 
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• the left hand side is non-negative and controls all spatial derivatives after 
applying the standard Hardy inequality (cf. 

• the second and third term on the right are essentially controlled by the 
E 2 [tjj] energy times the length of the time interval, (t 2 — Ti): 

g**** Vt* (QiV) Vt* ~ ^r 2 !^*^!} = < ex [d t ^] 

<C e -ei[flHfl+e-eaM (78) 

and e = 0. The fact that we need to borrow an e of e 2 [-0] is due to 
the fact that the t*<f> coordinates are not optimal near infinity. The cross- 
term g* $ term would decay much stronger in coordinates adapted to the 
asymptotically AdS end, which would allow us to estimate all terms by 
the weaker energy E 2 [ip]. 

• the first term on the right hand side is a boundary term, which can be 
estimated 



X>(ti,t 2 ) 



< sup / e 2 [ip] r 2 sm6drd6d(f) . (79) 



• the last term in (|77l) is controlled as previously by the last line in the 
energy estimate ([751) . 

It follows that integrating (|77|) furnishes the estimate 

r 2 sm6dt*drd9d<j)\r 2 \d r n t iP\ 2 + \^n^\ 2 



'V{T U T 2 )n{r>2B} 

< max(l, (r 2 - n)) 



e sup £ 2 [0] (t) + C* • £ 2 (0) 

re(n,T2) 



(80) 



Now note that 

e 2 [v] (<*) + ll^v-tll^o,-3 ( ^ n{r < 2R}) + iia^ii^ ds(St , n{r < 2H}) 

< C R ■ E 2 [ip] (t*) < C fl • £ 2 [0] (0) 



(81) 



follows right from the boundedness statement for the E 2 [ip] energy and estimat- 
ing the weights away from infinity. We can integrate (fSTj) in time and add it to 
([8U)) which yields (first without the boxed terms) 



E 2 m (t 2 ) + / E 2 [0] (r) dr < C M> i, a ,a • #2 M (n) 



■ max (1, T2 — ri) 



e sup £ 2 W(T)+C e -£ 2 M(0) 

t6(ti,t 2 ) 



(82) 
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The estimate also holds with the boxed terms included, as follows from adding 
(EBJ) and (HU). We claim that O implies E 2 [ip] (**) < E 2 [tp] (0) provided e is 
sufhciently small depending only on the parameters (the constant Cjif j a ) and 
leave the verification to the reader. 

Remark A.l. An easier proof is available if one is willing to go to H\ dS . The 
Carter operator 



commutes with the wave operator. Since &i and d 2 trivially commute, we have 



Ex [A S W'] (**) < E 1 [Q^j] (f) + Ex [dli] (<*) + Ex [%1>] (f) 

< Ex [QV] (0) + Ex [S%il>] (0) + Ex [d 2 t i] (0) (84) 



and we can control all derivatives on S 2 from controlling the Laplacian via ellip- 
tic estimates. This yields the desired gain, albeit at the level of three derivatives. 
This is analogous to commuting with angular momentum operators twice. 
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